Dissipation of the Rabi model beyond the rotating wave approximation: 
Quasi-degenerate qubit and uhra-strong couphng 
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Environmental influences on the dynamics of a coupled qubit-oscillator system are studied analyt- 
ically. Wc investigate the case of a quasi-degenerate qubit within the ultra-strong coupling regime 
for which the qubit frequency is much smaller than the frequency of the oscillator, and the coupling 
between the qubit and the oscillator is large, both of which invalidate the usually employed rotating 
wave approximation. In contrast to the standard quantum optics master equation, we explicitly 
take the qubit-oscillator coupling into account while microscopically deriving a dressed state master 
equation. Using the derived master equation, we discuss a spectroscopic technique which can be 
used to probe the dressed energy level structure of the qubit-oscillator system. 
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I. INTRODUCTION 



A two-level system (qubit) that interacts with a har- 
monic oscillator has been the cornerstone of various ex- 
perimental and theoretical studies [il,[2]- In many of these 
physical scenarios, the dynamics of the qubit-oscillator 
system is governed by the Rabi Hamiltonian Q : 



H 



Rabi 



haJo(Jz/2 + hoja^a + h/3uj{a + a^)ax. (1) 



In the historically most common experimental con- 
text of cavity QED, the coupling is extremely small (i.e., 
/? < 10~^) and the qubit and the oscillator are nearly 
resonant (i.e. wq ^ w) [l|. Under these conditions, 
it is valid to make an approximation of dropping the 
counter terms, acr_ and a'^a^ where ax = {(y+ + cr_), 
from the interaction Hamiltonian. This is the well-known 
Jaynes-Cummings domain of the rotating wave approxi- 
mation (RWA) , where remarkably simple exact solutions 
are available for eigen-energies and eigenvectors of HuaU 
[J| . The effect of the counter rotating terms become im- 
portant either because of ultra-strong coupling, |/3| > 0.1, 
or because of extremely large detuning, jcj — cjol ^ uj+ujo. 
Under these conditions, the RWA cannot be made 5]. 

Apart from the coherent dynamics generated by Hjiabi , 
there is incoherent evolution due to the system's unavoid- 
able interaction with its environment. The dissipative 
dynamics of the qubit-oscillator system has been inves- 
tigated extensively, focusing mostly on the regime where 
the RWA is valid [6l-[ll| . The incoherent evolution in this 
regime is phenomenologically modeled by the standard 
quantum optics master equation (SME) which assumes 
that the qubit-oscillator coupling plays no role in the dis- 
sipation mechanism. This approximation leads the SME 
to have non-physical predictions 12-16]. Although the 
predictions of the SME are not strictly correct, the errors 
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in its predictions are small when the coupling is small. 
However, the SME must be abandoned when the coupling 
is ultra-strong to invalidate the RWA. 

Developments in the area of circuit QED have made 
it possible to engineer systems that operate in regimes 
where the RWA breaks down |17H19| . Thus, the question 
of damping in the regime beyond the RWA has gained 
practical importance and has recently begun to be the- 
oretically explored [Ig, [2(]| - |23| . All these investigations 
focus on understanding the effects of increasing the cou- 
pling strength on the damped dynamics while keeping 
the qubit nearly resonant with the oscillator. 

Damping in the ultra-strong but far off-resonance 
regime has not yet been analyzed in the literature. The 
aim of this report is to explore this parameter space. 
In particular, we investigate the quasi-degenerate qubit 
regime where the qubit is far red-detuned from the os- 
cillator (wo <C Lo) and the coupling strength is ultra- 
strong [2j. As the SME becomes invalid in this param- 
eter domain, we adopt an alternate approach and fol- 
low Carmichael's treatment [IJI of investigating the dis- 
sipative dynamics of strongly coupled systems by micro- 
scopically deriving the composite system's dressed state 
master equation (DME). The DME approach has pre- 
viously been used to study the dissipative dynamics of 
ultra-strongly coupled qubit-oscillator systems while re- 
strictingattention only to the near resonance regime 
[la . |20| . |23| . Here we use the DME approach to explore 
the dissipative dynamics of ultra-strongly coupled qubit- 
oscillator system in the regime where the qubit is quasi- 
degenarate and is thus far away from resonance with the 
oscillator. 

In contrast to the SME, it will be shown that the 
DME predicts the correct steady state of the interact- 
ing qubit-oscillator system at thermal equilibrium. In 
the derivation of the DME, we model the qubit and the 
oscillator environments by separate bosonic baths. The 
effect of the energy level fluctuations of the qubit is also 
considered. Using the DME, we suggest a spectroscopic 
technique which can be used to probe the dressed en- 



ergy level structure of the Rabi Hamiltonian in the quasi- 
degenerate qubit and ultra-strong coupling regime. 



II. ADIABATIC APPROXIMATION 

In order to derive the dressed state master equation, 
we need to find the eigen-energies and eigenstates of the 
Rabi Hamiltonian in the parameter domain of interest 
where ujq <^ uj and the coupling is ultra-strong and thus 
the RWA becomes invalid. For this purpose, we use a 
form of an adiabatic approximation which was previously 
presented by Irish, et al. [23]. Here, we quickly sketch 
this approximation. 

Because wq is so small, the terms of zero-order and 
first-order importance in HuaU are completely different 
from those in the JC Hamiltonian. Here ^hajQCTz can be 
ignored in zero-order, and the remainder of the Hamilto- 
nian provides the eigen-energies and eigenstates forming 
the basis for the dynamical evolution. 

When ignoring the CTz term, the ax operator is con- 
stant, so we work in its eigenbasis and use m = ±1 to 
designate its eigenvalues. In its eigenbasis the zero-order 
Hamiltonian can be quickly diagonalized by factoring: 



Ho = hw{a* 

= nwb% - 



■ mP){a + m/3) — hjo0 , 



(2) 



where the new b operators are obviously displaced ver- 
sions of the a operators, and we have 



Ho\m,Nm) ^£N\m,N,ri), 



(3) 



where 

£n ^huj{N-p'^). (4) 

The eigenfunctions of Hq clearly have the form 

\m,N„,) = \m)(g)\N„,), 

where |?7i = ±) are the eigenstates of ax and \Nm) are 
oscillator states defined by: 

|iV™)=exp(-m/3(at_a))|iV), 

^D{-mP)\N), (5) 

where D{mP) is a displacement operator. The states, 
|A^), are number states of the oscillator: a^a\N) — N\N). 
We see from Eq. ^ that \N.m) are displaced Fock states, 
the displacement of which depends upon the state of the 
qubit. The energy eigenstates corresponding to the same 
N: | — , A^_) and |-|-, A^+), are degenerate. 

The CTz term in the Rabi Hamiltonian remains to be 
included. It promotes transitions among the Hg eigen- 
states, and it has non-zero matrix elements that need to 
be organized systematically. A straightforward approach 
is to write out the Rabi Hamiltonian matrix in the Hq 
eigenbasis, including all of the cr^ contributions, in order 



to show that almost all of them can be ignored on the 
basis of an adiabatic approximation that can be seen as 
a natural extension of the rotating wave approximation 
|25| . That is, fast terms are identified and eliminated 
in favor of slowly varying quasi-static terms that are re- 
tained. The full Hfiabi matrix, in stylized form, is dis- 
played here. Rows and columns of Hq are ordered in the 
sequence |-, 1_), |-h, 1+), |-, 2_), |-h, 2+), . . . , etc.: 



H 



Rabi 



£0 


GO 


— 


01 


— 


02 .. 


00 


fn 


01 


- 


02 


- . . 


- 


10 


£^ 


11 


- 


12 .. 


10 


- 


11 


£1 


12 


- . . 


— 


20 


— 


21 


£2 


22 .. 


20 


— 


21 


— 


22 


£2 .. 



(6) 



The diagonal elements are the bare energies given in (j4]) . 
As shown, there are two £n values for each N. The 
dashes in the matrix stand for elements that are zero, 
and the off-diagonal paired numbers are a simple short- 
hand for the matrix elements of the a^ term in the Rabi 
Hamiltonian. That is, 



MN=hojo{T,M^\a,\±,N±). 



(7) 



An adiabatic approximation is easily made by noting 
that every off-diagonal element in Hj^aU is associated 
with an oscillation at the frequency {M — N)u}. Only 
the diagonal elements and the off-diagonal terms labelled 
NN are static, i.e., are associated with a zero-frequency 
oscillation. In our adiabatic approximation, only these 
terms are kept, implicitly averaging all of the others to 
zero over a few periods 27r/w. That is, we discard all 
terms with time dependences that are assumed, in this 
adiabatic approach, to be too rapid to be perceptible - 
terms oscillating at the (high) oscillator frequency uj or 
any of its harmonics. 

The consequence of this adiabatic approximation is a 
reduced Hamiltonian in block-diagonal form, where the 
general iVth block is 2 x 2 dimensional. Each of these 
blocks is spanned by the states | — , NJ) and |-|-, A^+), and 
has the form: 



rriN) _( £n NN 



(8) 



where the off-diagonal elements can be evaluated exactly 
to (25J: 






(9) 



with Ljq{x) being Laguerre polynomials. The eigenstates 
and the corresponding eigenvalues of Hj^ are: 



*^r = 



.iri 



,N+)±\-,N' 



E^^£M±fi^{N4N+). 



(10) 



Earlier explanations and uses of the adiabatic approxi- 
mation can be found in |25l - l27t 

It should be noted that in determining the above eigen- 
values and eigen-functions, only the qubit frequency was 
assumed to be much smaller than the oscillator frequency 
and no approximation on the strength of the coupling, 
1/3 1, was made. In this report, we do not explore the 
entire parameter regime spanned by arbitrarily high cou- 
pling strengths but restrict our analysis to the regime 
where terms that are fourth or higher order in the cou- 
pling strength can be neglected, 0{f3'^) ^ 0. This focuses 
our analysis to the parameter domain which is experi- 
mentally feasible with current technol ogy or is likely to be 
realizable within the near future [13, [iMl ■ I^ this regime, 
one can approximate the exponential appearing in Eq. 
(ITTO to be: 



,-2,3 -^ 



l-2/?2 



(11) 



A modest upper limit on the coupling strength that sat- 
isfies the above approximation is |/3ma£c| =0.2. Coupling 
strengths greater than |/?| = 0.1 are considered to be 
ultra-strong 2J| as the validity of the RWA is known to 
break down under such strong interactions. Thus, the 
approximation, ©(Z?^) ~ 0, does not prevent us from 
exploring the ultra-strong parameter regime although it 
does restrict us from exploring the deep strong coupling 
regime where |/3| ^ 1 |24| . 

If we further restrict our analysis to the parame- 
ter regime where the oscillator excitation number obeys 
N <^ N,nax = 1/(2/3)^, the Laguerre polynomial appear- 
ing in Eq. ^ can be approximated to: 



LMi-iP^) w 1 - 4iV/32 



(12) 



Under the approximations given in Eqs. pil) and (|12p. 
the energy eigenvalues in Eq. ([10]) take the form: 



^N 






(13) 



il^iv)}' ^1^0 form a harmonic oscillator ladder with fre- 
quency aj_ and ground state energy — ?icjo(l — 2/3^)/2. 
We will refer to this oscillator as HO~ . 

In order to conveniently treat the two harmonic oscilla- 
tor structures, we define the annihilation operators, a±, 
and projection operators, 1±: 



i±= E i*^)(*^ 



4\ 



(15) 



In terms of these operators and within the adiabatic ap- 
proximation, the qubit-oscillator Hamiltonian becomes: 



Had =h(u+ala+ + ^(1 - 2/3^)1^ 

+ ft(cj_aLa_-y(l-2/32)l_ 



(16) 



From the above form of the Hamiltonian within the adi- 
abatic approximation, Had, it is clear that within the 
parameter regime for which loq -^ oj and 0{f3^) ~ 0, the 
qubit-oscillator composite system can be described by a 
set of two oscillators: "HO^ with frequency cj^ and HO^ 
with frequency w_ . 

We reiterate the important point that the validity of 
the adiabatic approximation relies upon the slowness of 
the qubit, ujq ^ uj. It imposes no restriction on the 
strength of the coupling parameter /3. However, the 
choice of analyzing only the parameter space for which 
O (/?■*) ~ 0, restricts our analysis to the regime where in- 
stead of using the adiabatic approximation, one could use 
the Schrieffer- Wolff (SW) transformation [2^] to find the 
eigen-structure of HjiaM ■ A comparison between the SW 
transformation and the adiabatic approximation is given 
in Appendix [Cj For the quasi-degenerate qubit regime, it 
is shown that the adiabatic approximation works better 
than the SW transformation. 



III. DISSIPATION 



where 



uj± = cli =F 2ujq(3 



(14) 



In the above equation, we have neglected the constant 
energy term, —hjjjS'^, appearing in £fq (see Eq. @). 



Note that ^+ 



E- = fia;o(l - 2/3^ - 4Af/32). Since 



the excitation number is restricted by A^ ^ 1/(2/3)^, for 
the parameter regime that we consider we always have 

We see from Eq. ([T3|) that for various values of N, the 
set of energies: {E^}, are equally spaced. This suggests 
that the corresponding set of eigenstates: (IvP^)}, form 
a harmonic oscillator ladder with frequency uj^ and with 
ground state energy hujQ{l — 2/3^ )/2. We will refer to 
this oscillator as HO^ . Similarly, the set of eigenstates: 



The dissipative dynamics of the qubit-oscillator sys- 
tem is usually studied using a standard quantum optics 
master equation (SME) [6-10]: 

p=-[H,p]+KV[a]p + jV[a^]p, (17) 

where p is the density matrix of the qubit-oscillator sys- 
tem and K and 7 are the relaxation rates of the oscillator 
and the qubit respectively that are coupled to separate 
zero temperature baths. The dissipator, V, defined as 

V[0]p^ hOpO'f ~0''0p-p0^0)/2, (18) 

generates non-unitary evolution of the system and the 
effect of the environment on the system dynamics are 



encoded in it. In the SME, the qubit-oscillator couphng 
parameter, j3, plays no role in the form of the dissipators 
and is introduced phenomenologicallym the Hamiltonian 
only to generate unitary evolution ]15|. Because of ne- 
glecting the coupling, the SME predicts phenomena that 
are physically incorrect J13 - ll6l |. For example, if the qubit 
and the oscillator environments are both at zero temper- 
ature, the Boltzmann distribution demands the steady 
state of the qubit-oscillator system be the ground state 
of the joint Hamiltonian, Hj^abi ■ In contradiction to this 
basic principle of thermodynamics, the SME incorrectly 
predicts that the system will evolve out of its ground 
state even when the baths are at zero temperature. 

In order to avoid the unphysical predictions of the 
SME, it is necessary to take the qubit-oscillator cou- 
pling into account while describing the damped dynam- 
ics. This is achieved by deriving the dressed state master 
equation (DME) of the interacting systems [IJI • 

In this section, we study the damping of the qubit- 
oscillator system operating in the quasi-degenerate qubit 
and ultra-strong coupling regime where the RWA breaks 
down. We provide a microscopic derivation of the dressed 
state master equation with the dressed states being given 
in Eq. pO|) . Note that these dressed states were de- 
rived using the adiabatic approximation and retaining 
the counter-rotating terms in the Rabi Hamiltonian. 

We model the reservoirs of the qubit and the oscillator 
by separate collections of oscillators with free Hamiltoni- 
ans: 



where 



H^sc = f^J2''^^l^^^ 



H, 



qubit 



-^E 



^J.\c\cx 



(19) 



The Hamiltonian governing the interaction of the qubit 
and the oscillator with their respective reservoirs is taken 
to be: 

A 

^,U.t = ^E'?^(4 + ^^)^- (20) 



The coupling strengths of the reservoirs with the oscilla- 
tor and the qubit are parameterized by hx and qx respec- 
tively. These coupling strengths are assumed to be small 
enough to allow us to make the rotating wave approx- 
imation and thus throw away the counter terms in the 
system-reservoir interaction Hamiltonian (see Eq. (|A7I) ). 
It also allows us to make the secular approximation and 
Born-Markov approximation when deriving the dressed 
state master equation in Appendix [A] 

When the reservoir is at zero temperature, the DME 
takes the form: 



p(t) = -^ [HAD,p{t)] + JoscP{t) + JgMtpit), (21) 



JoscPit) ^T{oj+)V[a+] p{t) 

-t-r(a;_)P[a_]p(t) 



N 



Jqutp{t) =^7((ijv)2? [I*w)(*^l] P{t)- 



(22) 



N 



In Eq. (|21|) , p{t) is the qubit-oscillator composite density 
matrix and Had is the Hamiltonian derived within the 
adiabatic approximation in Eq. (J16l) . The Liouvillians, 
Jose and Jqbit, in Eq. p2|) are super operators that en- 
code the effects of the reservoirs coupled to the oscillator 
and the qubit respectively. The damping rates associated 
with these reservoirs are r(a;) and 7(w). The expressions 
for these damping rates are given in Eqs. (jA16|) and 
(1X271) . 

We now try to understand the dynamical implications 
of the master equation ([2T|) . For doing so, we recall that 
the energy eigenstructure of the coupled qubit-oscillator 
system in this regime consists of two harmonic oscillators, 
T-LO , with respective frequencies: a;±. The energy eigen- 
states of each of these oscillators are dressed states in 
the sense that they are entangled qubit-oscillator states. 
Looking at Eq. ([^. we see that the super operators, 
Jose and Jqbit, consist of dissipators with three differ- 
ent types of lowering operators: a+, a_ and the set of 
operators {|^^)(^I'^|}- The term corresponding to the 
dissipator V [a+], induces damping within the oscillator 
HO . Similarly, the term corresponding to the dissipa- 
tor V [a-], induces damping within the oscillator 'HO~ . 
The respective damping rates of these oscillators are 
r(w-(-). The two oscillators are incoherently coupled by 
the term corresponding the the dissipator V [|^^)(^^|] 
which induces transitions from |^^) to |^^). The term 
V [|^E'jv)(^7vl] S^ts contribution from the reservoir con- 
nected to the oscillator and the qubit. Corresponding to 
the reservoir connected to the oscillator and the qubit, 
the damping rates for the transition |^^) — > l^/v) 
are T(u)i^) and ^(lom) respectively. These features are 
schematically shown in Fig. [T] 

When a system interacts with zero temperature baths, 
we expect on thermodynamic considerations that the sys- 
tem should equilibrate to its ground state. Thus, in the 
present scenario, we expect the steady state of the cou- 
pled qubit-oscillator system to be the minimum energy 
state 1^0 ). To find the steady state predicted by the 
DME, we have to find the state which does not evolve 
in time. Putting p{t) = \^q){^q\ in Eq. ([HI), we easily 
get that p{t) = 0. This implies that the steady state is 
\^q) and confirms our earlier statement that at thermal 
equilibrium, the DME correctly predicts the steady state 
to be the ground state of the coupled qubit-oscillator sys- 
tem. On the other hand, by putting p{t) = |^o )(^o | in 
the SME given in Eq. ([T7]) , it can easily be checked that 
p{t) 7^ 0. This incorrectly implies that the system will 
evolve even if it is initially in the ground state. 
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FIG. 1. Energy level diagram for the eigenstates, I'J'jv)) of 
the qubit-oscillator system in the quasi-degenerate qubit and 
ultra-strong coupling regime where the RWA breaks down. 
The dissipators, ©[a-] and ©[a+J, induce transitions within 
the oscillators "HO' and HO^ respectively. Transitions from 
HO^ to nO' is induced by the dissipator V [|*;^)(*+|]. 
In the long time limit, the steady state of the system is the 
ground state of the interacting qubit-oscillator system: 1*1/(7). 



As noted earlier, the discrepancy in the SME arises be- 
cause it assumes that the qubit-oscillator coupling plays 
no role in its dissipator. The dissipator of the SME pre- 
dicts the zero temperature steady state of the system 
to be \gO), where \g) is the ground state of the qubit, 
(cr^|e) = |e) and az\g) = — Iff))- We can quantify the 
difference between the predictions of the SME and DME 
by calculating the distance, d, between the states I'J'q") 
and \gO): 



d=l-\{%\gO)\, 



(23) 



We see from the above equation that as the coupling 
strength becomes small, the distance d tends to zero and 
thus the prediction of the SME starts to coincide with 
the prediction of the DME. This agrees with the com- 
mon notion |12h16| that the applicability of the SME is 
restricted only to the parameter domain where the cou- 
pling between the qubit and the oscillator is small. 

Similar to the derivation of the master equation in the 
zero temperature case, the master equation for the case 
when the reservoirs are at finite temperatures can also 
derived. For simplicity, we do not explicitly discuss the 
finite temperature case here and only give the results. If 
the qubit and the oscillator baths are at temperature T, 



the steady state of the DME can be calculated to be 



Pss = 



Xr Ig-HAn/kBTl ' 



(24) 



where fc^ is the Boltzmann constant. This steady state 
correctly predicts the occupation probabilities of the 
dressed states of the coupled system to be the canonical 
Boltzmann distribution. It can also be shown that the 
finite temperature DME satisfies the principle of detailed 
balance whereas the SME does not. 



IV. QUBIT ENERGY LEVEL FLUCTUATIONS 

The energy level splitting of a two level system can be 
sensitive to the environment with which it interacts. The 
effect of energy level fluctuations can be accounted for by 
introducing the following fluctuating Hamiltonian: 



Hfit)^hf{t)a, 



(25) 



where f{t) is a stochastic term. Thus, in the presence of 
energy level fluctuations of the qubit, the total Hamilto- 
nian governing the dynamics of the qubit-oscillator sys- 
tem becomes: 



H{t) = H^ 



AD 



Hfit). 



(26) 



With the above Hamiltonian, it is now a question as 
to what is the ensemble averaged effect of Hf {t) on the 
evolution of the coupled qubit-oscillator system. In Ap- 
pendix [Bl we find that in the presence of the fluctuating 
Hamiltonian, the evolution of the system is governed by: 



m 



ih 



[Had, pit)] 



If 



[5,,[^„p(i)]], (27) 



where we have converted back to the Schrodinger picture 
from the interaction picture that was used in Eq. (JBIII) . 
The operator, S'^, is the dephasing operator which has 
the form: 

S. = Y,{N+\N^) {\^%){^%\ - |*^)(*^|) . (28) 

N 

To understand the dynamical implications of Eq. (|27|) , 
we look at the evolution of the various matrix elements 
of the coupled qubit-oscillator density matrix: 



(*^|p(t)l*M) 



X (*±|p(i)|*± ), 



{M^\M-)Y 



= --^0(/34)(*±|p(t)|*±), 

{n\KmM) = - Y {{N+\N-) + {M+\M.)f 

X {n\pit)\Hi)- 



(29) 



(30) 



We see from Eq. ((291) that the population of the vari- 
ous energy eigenstates of the qubit-oscillator system does 
not change because of the energy-level fluctuations of the 
qubit. This is in contrast to what is found in the analysis 
of qubit-oscillator evolution in the near resonant ultra- 
strong coupling regime where energy-level fluctuations of 
the qubit tend to induce population transfer between the 
eigenstates of the Rabi Hamiltonian 16]. We also notice 
from Eq. (|29l) that qubit level fluctuations does not have 
any dynamical consequence if the population is entirely in 
the harmonic oscillator HO^ or entirely in 'HO~ . From 
Eq. ([50)) we see that any initial coherence between the 
various eigenstates oi'HO'^ and "HO", i.e. between states 
\^%) and \^m)^ S^t^ exponentially damped to zero. 



V. DRIVEN DAMPED QUBIT-OSCILLATOR 

SYSTEM 



We now consider the case when the oscillator is 
pumped by a classical drive of frequency ujp and ampli- 
tude Qp. The Hamiltonian modeling the pump is: 

Hp = hQpia + at)(e*"p* + e"*""*), 

= hnp{S + S'')ie"^''' + e~"^'-*), (31) 

where S is the dressed lowering operator defined as 



k>] 



a_+a+-2/3^|vI/^)(vI/+| 



(32) 



N 



In the above equation, wc have used the matrix elements 
of [a + a^) derived in Eq. (|A21|) . 

Including the pump Hamiltonian and the effect of qubit 
energy level fluctuations in Eq. (pTj). we get the master 
equation describing the evolution of the pumped qubit- 
oscillator system to be: 

p{t) =1 [HAD,p{t)] + ^[S + S\p{t)\ ie^-^' + e-^"-*) 
+ JoscPit) + Jgutpit) - y [^., [S.,pm , (33) 

In order to remove the explicit time dependence from the 
right hand side of Eq. ([55]) , we go to a frame rotating at 
the pump frequency and define the density matrix in the 
rotating frame by x(i)- 



x{t) 



iujp{a:a-\-~\-a_a-)t 



p{t)e 



'iujp{a,a^-\'a_ a—)t 



(34) 



In this rotating frame, the term ^ |'I'^)('I'^| in the driv- 
ing part of the Hamiltonian rotates at the pump fre- 
quency, LOp. In the spirit of the rotating wave approx- 
imation, we average this term to zero. Two other terms: 



Se '"^p* and S'^e^'^p*', also average to zero. Thus the mas- 
ter equation in the rotating frame becomes: 



N 



If 



[s,,[Sz,xm 



(35) 



where 



+ TV[a±]x{t), 

A± =iO± — LOp. 






a± -f a^,x(i) 



(36) 



In the above equations, we have assumed that r(ci;-)_) — 
r(a;_) = r and ^{ujn) + '^0^^{^n) = k does not depend 
on A^ for simplicity. The Liouvillians C^ and C~ gener- 
ate evolution of the harmonic oscillators TiO^ and TiO^ 
respectively. The second-last term in Eq. ([55]) transfers 
population incoherently from TiO^ to TiO^ at the rate 

The structures of the Liovillians, £+ and £^, are for- 
mally the same as that of a Liouvillian that generates 
evolution of a bare driven damped harmonic oscillator 
[23]. The steady states generated by C^ are coherent 
states, |a±), and are easily evaluated to be: 



l«±) 



\N 



Y^ ("^ \^± 



N 



'N\ 



Nh 



i L^ 



a± = 



(37) 



The fact that \a±) are indeed the steady states of £^ can 
be checked by verifying that C^\a±) = 0. We reiterate 
the important point that these steady states are coherent 
states in the dressed basis, {|*^)}. 

The long term behavior of the system dynamics can be 
understood by looking at the evolution of the following 
terms: 



N N 



N 



N 



In the above expressions, we have used Eq. (|29|) . The left 
hand side of Eq. (|38|) corresponds to the total population 
of the system within the "HO^ oscillator. Similarly, the 
left hand side of Eq. (pg]) corresponds to the total pop- 
ulation within the T-LO^ oscillator subspace. From Eqs. 
we see that in the long time limit, 



hm Y.^^N\x{mi) = 0. 

N 



tK» 



(40) 



Thus, in the steady state, there is no population in the 
osciUator, "HC*". From Eq. ([M]) . we see that ah the 
initial population in T-LO is incoherently transferred to 
HO~ . We conclude from these considerations that in 
the limit, Ki ^ 1, all the population is found only within 
the HO^ oscillator subspace. As noted in Eq. ((37|). the 
steady state within the HO~ oscillator subspace is |a_). 
Thus in the long time limit, any initial state will evolve 
to the steady state 



Xss = |a_)(a_|. 



(41) 



Defining the sum of average excitation numbers of the 
HO' and HO^ oscillators to be 

Af{t) ^ {aUt)a-{t)) + {al{t)a^{t)), 

= N-{t)+N+{t), (42) 

we see from Eq. (pij) that in the steady state, we have 



^fss^^f:. 



VI. SPECTROSCOPY 



(43) 



In this section, we will see how one can use the sum 
of average excitation numbers, Mss, to experimentally 
find out the energy eigenstates of the composite qubit- 
oscillator system. 

The bath modes coupled to the system allow the 
system energy eigenstates to be probed through spec- 
troscopy. In the physical scenario where the oscillator 
bath mode is used as a measurement channel, the out- 
put signal will depend upon the oscillator bath opera- 
tor: b\{t). From the system-bath interaction Hamilto- 
nian given in Eq. (jA7p . it is clear that the oscillator bath 
operator evolves as: 



bx{t)e'''-* =6a(0) 



ihx I dt' S{t)e 
/o 



'^c*'^A*' 



(44) 



where S{t) is defined in Eq. ((32)) . 

The rate at which the system dissipates energy into the 
oscillator bath is an experimentally measurable quantity. 
If the bath is initially in the vacuum state, the power 
that is dissipated by the system into the oscillator's bath 
in the steady state is given by 

Pss = lim , 7 E ^i^A(fcl(t)&A(i)).., 

p- /'OO 

= , l^JJLi 7 / dvg{v)h^{v)v 

X f dh f dt2{S\h)S{t2))sse-'''^'--''\ (45) 



The sum over modes, A, has been converted into an inte- 
gral with giy) being the density of states. Let us suppose 
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FIG. 2. Two-tone spectroscopy to probe the |^^) O |^jv) 
transitions. The dashed (green) arrows indicate the pump 
frequency which is taken to be resonant with HO" , uip = u)-. 
This pump is off-resonant with the T-LO'^ oscillator, cjp 7^ lo+. 
The solid (blue) arrows indicate the spectroscopy drive which 
is scanned around uq. Whenever ujs ~ liiv, population gets 
transferred from |^^) to IvP^). 



that the oscillator is driven by a classical pump of fre- 
quency ijp which is close to the bare oscillator frequency 
id. In this case, as in the previous section, the term 
^|^^)(^^| in the operator S{t) oscillates at Wp and 
can be averaged to zero. Following the usual approach 
[30J of using the Markov approximation to calculate the 
integral in Eq. ([^ . we get: 

Pss = Swp — -^ f (a^a+)ss + (aLa_)ss j , 
T{ujp) 



huJr, 



-K, 



(46) 



From the above equation, we see that the experimentally 
measurable quantity, P^s, is proportional to the sum of 
the steady state average excitations of the 'HO~^ and the 
HO^ oscillators. 

We now discuss how one can use Pgs to probe the en- 
ergy spacing between the dressed states |^^) and |^^), 
i.e. states having the same values of N. For this purpose, 
we propose a spectroscopic technique that is similar to 
the one discussed in 31]. The technique requires driving 
the oscillator by a pump drive with frequency ujp and a 
spectroscopy drive with frequency Ws . The pump drive is 
taken to be resonant with the oscillator HO^ , ujp = lo-, 
and the frequency of the spectroscopy drive is scanned 
near the qubit frequency, uJs ~ wq. The setup for this 
two-tone spectroscopy is shown in Fig. ^. 

In the absence of the spectroscopy drive, the pump 
drive forces the system to oscillate and in the steady 
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FIG. 3. Numerically evaluated percentage reduction in the 
steady state excitation number, J\fas, as a function of the spec- 
troscopy drive frequency, ujs- The pump frequency is taken to 
be resonant with the "HO^ oscillator, Up = uj-. The qubit fre- 
quency and the coupling parameter is taken to be ojq = O.Soj 
and j3 = 0.1 respectively. The damping rates are taken to be 
r = 6k = 37/ — 3cJo/? /5. The spectroscopy drive strength 
is taken to be Q.s = k. Maximum reduction in the power is 
achieved whenever uq ~ u>s = 2ujo/3'^{2N + 1). The various 
curves correspond to various pump strengths which are pa- 
rameterized by the corresponding average excitation numbers, 
Af^a, calculated in the absence of the spectroscopy drive. For 
higher average excitation numbers, peaks with higher values 
of A'^ in Eq. H48p start getting more dominant (this is easily 
seen in the figure for (ujq —u>s)/2ljoI3'^ ~ 5), whereas peaks at 
lower values of A'' get smaller. 



state, when all the population relaxes to the HO^ os- 
cillator, the average excitation number is given by Eq. 
((33]) with A_ = (since ojp = a;_): 



M,s ^ M-, ^ AQl/r^ 



Now in the presence of the spectroscopy drive, when- 
ever the |^]v) ^^ l^jv) transition becomes resonant with 
the spectroscopy drive frequency, i.e. when 

Lj, = c^o(l - 2/3^ - 4iV/32), 



2l^o/32 



= 2iV + l, 



(48) 



appreciable population gets transferred from the HO' 
oscillator to the HO'^ oscillator. This causes a reduc- 
tion in the excitation number of T-LO^ . Since the pump 
drive is off-resonant with the T-LO^ oscillator, Up ^ a;+, 
the population that gets transferred does not get reso- 
nantly driven. Hence, the reduction in the excitation 
number of the T-LO~ oscillator does not get compensated 
by the increase in excitation number of the HO oscilla- 
tor. This results in a net reduction of Afss whenever the 
spectroscopy drive frequency, w^, satisfies the resonance 
condition given in Eq. (|48l) . Reduction in the average ex- 
citation number manifests itself through a reduction in 
the dissipated power. Thus, by monitoring the decrease 
in the dissipated power, Pgs, as the spectroscopy drive 
is scanned around ujq, one can determine the various fre- 
quencies, Wat = a;o(l — 2/?^ — AN(3^), that are associated 



with the transitions between the dressed states |^^) and 

In order to include the effect of the spectroscopy drive, 
we add to the master equation (j33p a spectroscopy drive 
Hamiltonian, H^, which is identical to the pump Hamil- 
tonian, Hp given in Eq. pip , but with flp and ujp re- 
placed by spectroscopy drive strength, fig, and frequency, 
LOs, respectively. We then numerically solve the master 
equation and evaluate the steady state average excitation 
number Afss- We plot in Fig. |3]the percentage reduction 
in J\fss as a function of the spectroscopy drive frequency, 
Ldg- In agreement with our understanding, we clearly see 
from the figure that whenever the resonance condition 
defined in Eq. (|48l) is satisfied, there is a reduction in 
the the average excitation number. This provides a clear 
signature using which one can probe the transition ener- 
gies between the states |^^) and |^E'^). 

In Fig. [3l we also study the effect of the strength 
of the pump drive on the reduction of the steady state 
power. Using Eq. (|47l) . we see that the pump strength, 
flp, can be parameterized by M^^, which is the average 
excitation number of the HO^ oscillator in the absence 
of the spectroscopy drive. Since the pump is taken to be 
a coherent drive, one expects that in the steady state, 
the populations in the |^^) states will have a Poisson 
distribution. As Af~g increases, population in the higher 
I'^Jf) states also increases. This increases the ability of 
the spectroscopy drive to transfer population to the |^^) 
states with higher values of TV. Thus, we expect the peaks 
with higher values of N to be more dominant as N~ 
increases. This is in agreement with the plots in Fig. [3] 
where we see that as Af~ increases, the peaks with higher 
('4y^ values of N in the transition |^^) — ^ I^tv) start getting 
more pronounced while peaks with lower values of N get 
smaller. 



VII. REGION OF VALIDITY 

The eigenvalues, eigenfunctions and the Lindblad mas- 
ter equation derived in this report are based on certain 
approximations. For a better understanding of the valid- 
ity of these analytic formulas, we list below the various 
constraints under which these results were derived: 

(a) u!q 'ti Lo: Quantifies the quasi-degeneracy of 
the qubit. Under this condition, the time scales 
associated with the qubit evolution and the oscil- 
lator evolution are well separated and one can con- 
clude that the fast oscillator adiabatically follows 
the slow qubit and the fast oscillations of the qubit 
at or above the oscillator frequency can be aver- 
aged to zero. This validates the use of the adia- 
batic approximation (Sect. |TT|. The maximum pos- 
sible qubit frequency within which the adiabatic 
approximation is valid is ojq = O.'Soj [27|. 

(b) 0{I3^) ^ 0: Justifies the approximation made 
in Eq. (|lip . Even within this approximation, the 



coupling can be strong enough, |/3„ 
invalidate the RWA. 



0.2, to 



(c) N ^ (1/2/3)^: Allows us to restrict the power 
series expansion of Ln to the first two terms only. 
It is used in Eq. P^ and sets a limit on the validity 
of relations ([T5|) . 

(d) Secular approximation: Assumes that the en- 
ergy difference between dressed states is much big- 
ger than the relaxation rates between them (see 
Eq. (jAlSj). We have also used this approximation 
to remove the "counter rotating terms", B{t)S{t) 
and B\t)S''{t), from the system-bath interaction 
Hamiltonian to arrive at Eq. (|A7|) . 



(e) Born approximation: Assumes that the cou- 
pling of the qubit and the oscillator to their re- 
spective baths, qx and hx, are weak. This allows 
restricting the system-bath interaction to second 
order in the coupling strength. Also within this ap- 
proximation, one assumes that the reservoir state is 
not effected due to its interaction with the system. 
Born approximation is used in Eqs. 



(f) Markov approximation: Assumes that the bath 
correlation time is much smaller than the time scale 
in which the system evolves. Markov approxima- 
tion is used in Eqs. (jAlip and (jAlsp . 



VIII. CONCLUSION 



quasi-degenerate qubit regime. It was shown that the 
eigenstructure of the qubit-oscillator system comprises 
of two sets of harmonic oscillators, HO^ . These oscil- 
lators have different frequencies, lu± = oj ^ 2ajo/3^, and 
different ground state energies. 

The damping rates, T{uj±), T{ujn) and j{ujn), enter- 
ing the DME depend upon the bath density of states 
evaluated at the eigen-frequencies of the coupled qubit- 
oscillator system. Since the eigen-frequencies depend 
upon the coupling strength, /?, it is in principle possi- 
ble to change the damping rates just by tuning the cou- 
pling. This dependence of damping rates on the coupling 
strength does not arise within the SME. 

The master equation governing the dynamics of the 
damped qubit-oscillator system when the oscillator is 
driven coherently is derived in Sec. |Vl When the drive 
frequency is taken to be nearly the same as the bare os- 
cillator frequency, the steady state of the driven damped 
system is found to be a coherent state in the subspace 
spanned by the oscillator, HO^ . 

A two tone spectroscopic technique which can probe 
the energy differences between the eigenstates |^^) and 
I^P^) is discussed in Sec. I VII By solving the dressed mas- 
ter equation, we find that whenever the spectroscopic 
drive frequency becomes resonant with the I'l'^) — > 
IvP^) transition, the steady state power dissipated by the 
system, P^s, reduces. This provides a tool with which 
one can probe the eigenenergies of the composite qubit- 
oscillator system in the quasi-degenerate qubit and ultra- 
strong coupling regime where the RWA breaks down. 



In this report, we investigated the effect of the envi- 
ronment on the dynamics of a coupled qubit-oscillator 
system. We restricted our analysis to the parameter 
regime where the qubit is quasi-degenerate, wq ^ uj, 
and the coupling is allowed to be ultra-strong. For such 
a parameter regime, the rotating wave approximation 
breaks down. Also in this parameter domain of interest, 
the standard master equation which is commonly used 
to describe the effect of the environment on the qubit- 
oscillator system is known to fail. For example, when 
the baths are taken to be at zero temperature, the SME 
incorrectly predicts the ground state of the composite 
qubit-oscillator system to be not the equilibrium state. 

To get around the nonphysical predictions of the SME, 
we microscopically derived the dressed state master equa- 
tion. The DME is shown to predict the correct thermal 
equilibrium state. We showed that the differences be- 
tween the predictions of the SME and the DME decreases 
as the coupling strength decreases. This is in agreement 
with the fact that the SME is valid only when the cou- 
pling is small but must be abandoned when the coupling 
is ultra-strong. 

For deriving the DME, we needed to find the eigen- 
energies and eigen-functions of the Rabi Hamiltonian. 
We achieved this by using the adiabatic approximation. 
The adiabatic approximation is valid only within the 
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Appendix A: Master equation 

In this appendix, we will derive the master equation for 
the qubit-oscillator system when the qubit and the oscil- 
lator are interacting with separate bosonic baths. The 
analysis will be based on Born-Markov and secular ap- 
proximations and will be carried out in the dressed basis. 
We will assume that the initial state of each bath is a 
thermal state. We explicitly carry the calculations only 
for the case when the bath is at zero temperature. 

We first consider the case in which only the oscillator 
is coupled to its corresponding bath. For this setup, the 
system-bath interaction Hamiltonian is given by: 



HL^hJ2hxibl + bxn 



(Al) 



where C is a system operator that, according to Eq. (|20p . 
is (a^ -t- a) . We now write the interaction Hamiltonian in 
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the interaction picture: 



HL{t)^fi E /iA(6ie-^* + feAe-''^^*)Qfe|j)(fc|e-^'=^*, 



A,b-,fc) 



(A2) 



where the set of states {\j)} are the eigenstates of the 
coupled qubit-oscillator system as given in Eq. ([T(I| . The 
coefficients, Cjk, and the frequencies, A^j, are defined as: 



In Eq. (|A8[) . p{t) is the reduced density matrix of the 
qubit-oscillator system in the coupled system interaction 
picture, Tr^ refers to taking the trace over the bath de- 
grees of freedom and ps(0) is the initial density matrix 
of the bath. For the zero temperature result, pb{0) is 
taken to be the vacuum state. Putting the interaction 
Hamiltonian, Eq. (jA7|) . in Eq. (|A8P and collecting all 
the non zero terms that contribute to the evolution, we 
get: 



Afej = ujk ~ ojj 
Cjk = {j\C\k). 



(A3) 



m 



1 



m^ Jo 



dt'[S\t)S{t')p{t'){B{t)B\t')) 



In order to aid making the secular approximation, we 
write the interaction Hamiltonian in a more convenient 
form: 

X,(k>j) 



-ivxt 



)C,u\3){kV 



-iAkjt 



+ h Yl hxible'^-' + bxe- 

\,(k<j) 

+ fi E h.ible^'^^' + 6,e— ^*)Q,b-)(j|. (A4) 

Taking C — {a + a''), it can be easily shown that Cjj — 0. 
So, we get: 



+ H.C, 

^iBit)+B^it))iSit)+SHt)), 
where we have defined 



B{t) = hJ2hxbxe-'''-', 

A 

{k>j) 



-iAkit 



(A5) 



I A kit 



(A6) 



It should be noted that B{t) is the lowering operator for 
the bath and S{t) is the analogue of the system lowering 
operator in the dressed basis. According to the first sec- 
ular approximation, we neglect the terms B{t)S{t) and 
B{t)''S{ty from the interaction Hamiltonian to get: 



HUt)^Bit)SHt) + B\t)S{t). 



(A7) 



Within the usual Born approximation of treating the 
system-bath interaction only up to second order in the 
coupling strength, hi, and assuming the reservoir state 
to be unchanged during the dynamics, the evolution of 
the reduced density matrix of the qubit-oscillator system 
in the interaction picture follows the following integro- 
differential equation: 



m 



Trg 

(ih)^ 



dt' 



Hi(t) Hj{t'),p{t')<E>pBiO) 



(A8) 



- Sit)pit')SHt'){B{t')B^{t)) + h.c}j , (A9) 

where {B{t)B^t')) = Tib {B{t)B^t')pB{0)} is the bath 
correlation function. From Eq. (|A6|) it can easily be 
seen that {B{t)B^t')) = {B{t - t')B^O)). Let us now 
evaluate the first integral: 

I,{t) ^^ j' dt'S\t)S{t')p{t'){B{t)B\t')), 

^ f dTS\t)s{t - Typ{t t){b{t)b\q)). 



h^ 



(AlO) 



Based on the assumption that the bath correlation 
function, (_B(r)_B'l'(0)), is peaked only for an interval 
T — tb which is much smaller than the relaxation time of 
the system, r^j, one can make the Markov approximation 
and replace p{t — r) in Eq. ()A10p by p{t) to get: 

^iW ^^f dTS\t)S{t~T)p{t){B{T)B\Q)). (All) 



Using Eq. (jA6|) . we get: 

h{t)= E C,kCU^)m{k\p{t)fokUt), (A12) 



where 



^ JQ 



i{ux-Akj)T 



(A13) 



Converting the sum over the various bath modes in the 
above expression to an integral, we get: 



Jjkhnit) = e^(^"'-^''^)* / dvg{v)h\v) 
Jo 

X f dre-'^"-^"^^'' , (A14) 

where g{i') is the density of modes for the bath. Assum- 
ing the relaxation time of the system to be much longer 
than the correlation time of the bath, the upper limit of 
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the T integral can be set to infinity. Neglecting the di- 
vergent frequency shift, fjkj'k'{t) can then be evaluated 
to get: 



fjklmit) = £%)e'(A„.~A.,)*, 



where we have defined 



r(Afe,) = 27:g{Ak,)h^{Akj) 



(A15) 



(A16) 



the usual Fermi Golden Rule result. 

The terms in Eq. ()A15|) that are oscillatory are ad- 
ditional non-secular terms and in the second secular ap- 
proximation, these terms are averaged to zero. Thus any 
term, fjkimit), for which Ami 7^ Akj is neglected. The 
secular approximation implies: 

fjkim = ^^^5{Ami - Akj), (A17) 

where 6{Am.i — Akj) stands for the Kronecker delta func- 
tion. This approximation is justified only if the relaxation 
rate of the system is much smaller than the oscillatory 
frequency ;22i, ^], i.e.: 



r(Afcj)< \Aml-Akj\. 



(A18) 



Thus, under the Born-Markov and secular approxima- 
tions, one arrives at the following form for Ii(t): 



hit)= J2 C,kCUm){l\j){k\p{t) 



(A19) 



We now look at the explicit form of the matrix ele- 
ments, Cjk- In the quasi-degenerate qubit regime, the 
eigenstates of the qubit-oscillator system are: 



I* 



L(i 



For C ~ a + a\ we get: 



N.)±\~,N' 



(A20) 



-M/' 



= VmSn,ai~i + VM + ISn,m+i, 



(A21) 



Putting Eq. (IMTj) in Eq. ((XT9)) . we get: 

Ii{t) — — I r(w_|_)a^a_|_ + r(a;_)a_a_ 



4/32^r((iA,)|vI/+)(vI/+|)p(t). (A22) 



N 



where wat = wo(l - '^NjS'^ - 2/3^). 



Following the same arguments used to evaluate the in- 
tegral /i(i), one can now evaluate the rest of the integrals 
appearing in Eq. ()A9p . In particular we have: 



hit) 



ft2 



dt'S{t)p{t')S\t'){B{t')B\t)), 



= 2(r('^+)a+PWa+ +r(w-)a-/5(i)aL 

+ 4/?2^r(^A.)|*^>(*+|p(i)|vI/+)(v|/^|). (A23) 

N 

Putting these integrals in Eq. (jA9l) and going from the 
interaction picture to the Schrodinger picture, we get: 

p{t) ^ ^[H,p{t)]+V{uj+)V[a+]p{t) 
in 

-|-r(w_)P[a_]p(i) 



N 



ih 



[H,p{t)]+JoscP{t): 



(A24) 



where the dissipator V is defined in Eq. ()T8|) . The Lind- 
bladian, Josa describes what effect of the oscillator being 
coupled to a zero temperature bath. 

In order to derive the evolution equation for the system 
when the qubit is also coupled to a bath, we note that 
the matrix elements of ax in the energy eigenbasis are: 



= 0, 

Cn_m_ = Cn+Mj^, 
Cn+m^ = 5nm, 



(A25) 



Using Eq. (|A25|) and under the Born-Markov and secular 
approximations, the master equation now takes the form: 



p{t)^-T:[H,p{t)]+JoscP{t) 
in 



J2^{C.^)V[\^-^){^+\]p{t), 



N 

- [H, p{t)] + JoscP{t) + JqMtp{t), (A26) 



In the above equation, the damping rate 7(0;) is defined 
as: 



7(0;) = 2Tig'{uj)q^{uj), 



(A27) 



where g'{(jj) is the density of states and q{(jj) is the cou- 
pling strength of the qubit with the bath evaluated at 
the frequency lo. The master equation (|A26P describes 
the system evolution for the case of the qubit and the 
oscillator coupled to zero temperature baths. 
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Appendix B: Energy level fluctuation of qubit 

To understand the effects of fluctuating energy-level 
splitting of the qubit, we look at the evolution of the 
qubit-oscillator density matrix generated by the Hamil- 
tonian: 



H{t)^HAD + Hf{t). 



(Bl) 



Going to the interaction picture, the evolution equation 
has the form: 



Pfit) = ^[^^{t),pf{t)], 



(B2) 



where pf{t) = e""^"* pf{t)e'~'"''"^ and ^^(i) = 
^iHAot^^g-iHAot QJ.Q operators in the interaction picture. 
The subscript / in Pf{t) indicates the evolution of the 
density matrix for a particular realization of the fluctu- 
ating term f(t). 

Formally integrating Eq. (|B2p . we get the following 
integro-differential equation: 

P/(i)=-*/W[^.(t),P/(i)] 

'dt'f{t)f{t')[a4t),[a4t'),pf{t')]]. (B3) 

Assuming that the fluctuations are small, we only con- 
sider the effect of f{t) up to its second order. Moreover, 
we assume that the fluctuations of f{t) is not effected 
by the dynamis of the qubit-oscillator system. These are 
equivalent to the Born approximations. Taking the en- 
semble average of the various realizations of f{t), we get: 



m =pf{t), 



*/W ^zit),Pfit) 



- / dt7(i)/(i') ^.(i), ^z{t'),pf{t') 

Jo L L 



(B4) 



where the overhead bars indicates ensemble averages. 

The statistics of f{t) is in general governed by the 
physical system under consideration. For simplicity, we 
assume that /(i) has zero mean. Moreover, we assume 
that it is a stationary random process having the fol- 
lowing spectral decomposition of its two time correlation 
function: 

{fit)fit')) = {fit-t')fm, 

= — / diygiiy)e~"'<^'-*\ (B5) 

With these considerations, the evolution equation be- 
comes: 

-1 j-t 1-00 

m^-TT- / dt' / dz^^He— (*-*') 

x[d.{t),[a.{t'),p{f)\] (B6) 



We now find out the matrix elements of dzit) in the 
eigenbasis, |^^), of the adiabatic Hamiltonian, Had- We 
find that 



{^%\dM\Kl)N^M 

\^s)\n) 



(*^i" 






(*±|a.(i)|*±)=.±(iV+|7V_ 



(B7) 

(B8) 
(B9) 



In the evolution Eq. (jB6p . the operator az{t) appears 
quadratically. This means that the contribution of terms 
in Eqs. ((B7)) and (|B8)) wfll be of fourth order in the 
coupling: O (/3'*) , and hence we set them to zero. With 
this in mind, we can approximate dz{t) in Eq. (|B6P to 
be Sz, where: 

Sz = E(^+l^-) (l*^)(*^l - l*iv)(*^l) ■ (BIO) 

N 

It should be noted that Sz is time independent. Under 
the Markov approximation, the evolution Eq. (jB6|) can 
now be evaluated to be: 



m^-^-^[s.,[s.rpm^ 



(Bll) 



where 7/ is the dephasing rate quantified by the zero fre- 
quency component of the two time correlation's spectral 
distribution: 



7/ - 5(0). 



(B12) 



Appendix C: Adiabatic approximation and 
dispersive regime beyond RWA 



As mentioned in Section FVlIl our analysis is restricted 
to the quasi-degenerate qubit regime, wq < O.Sw, and 
where the coupling is no stronger than |/3ma2;| =0.2. This 
implies that all our analysis is valid only in the range 
where the coupling is much smaller than the detuning: 
a; 1/3 1 <C I a; — Wo I and the detuning is comparable to the 
bare qubit and the oscillator frequencies: io — loq r^ \lo + 
ujo I . This parameter regime is called the dispersive regime 
beyond the RWA [H. 

It was shown by Zueco, et al. [23, that the eigenvalues 
of Hiiabi in this parameter regime can be approximately 
found by using the Schrieffer- Wolff (SW) transformation. 
According to the SW transformation, one unitarily ro- 
tates Hjiabi to an approximate diagonal Hamiltonian: 
Hsw = G~^ HfiabiG^ , whcrc S is an anti-hermitian ma- 
trix (see Eq. (13) of ^]), to get: 

Hsw —^0-^ + hcua^a + hujo{^ -](Tz(a^ + a)^. 

2 Vwq — lo^ ' 

(CI) 

Note that unlike in the dispersive regime within the RWA 
[2], the oscillator excitation number, o^a, is not a con- 
stant of motion. Corresponding to the two eigenvalues of 
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FIG. 4. Fidelity of the first A'' eigenfunctions calculated 
within the adiabatic approximation (red crosses X) and the 
Schrieffler- Wolff transformation (black bold dots •). When 
wq <C w, the adiabatic approximation has better fidelity. 
When Wo ^ w, the adiabatic approximation breaks 
down. 



az , the eigen-energies of Hsw can easily be calculated to 
be 



^N 



±fia;o/2 + feI'±(iV + l/2), 



where 



UJ± = cj I 1 =p 



4cjwo^2 s 1/2 



(C2) 



(C3) 



Unlike the adiabatic approximation, which works only 
if the qubit frequency is much smaller than the oscillator 
frequency, the validity of the SW transformation does not 
depend up on the sign of the detuning: wq — w. Thus, 
Eqs. (jC2p and (|C3I) are valid for a red detuned qubit, 
Wo < w, as well as for a blue detuned qubit, uiq > uj. In 
the quasi-degenerate regime, we have wq ^ w and so we 
can approximate the frequencies uj± in Eq. (jC3p to: 



uj± — uj^ 2uJof3 



(C4) 



Using the above approximate expression for Cj± in Eq. 
(|C2p . and comparing it with Eq. ([TU]) . we see that the 
eigen-energies calculated using the adiabatic approxima- 
tion matches exactly with the eigen-energies derived us- 
ing the SW transformation: E^ — E^. 

Although the eigen-energies predicted by the adiabatic 
approximation and the SW transformation are the same, 
we now show that the eigenfunctions calculated within 
the adiabatic approximation have better fidelity with the 
exact eigenfunctions than the ones calculated according 
to the SW transformation. The eigenfunctions calcu- 
lated within the adiabatic approximation are given in 
Eq. (|10|) . For the SW transformation, the eigenfunc- 
tions are numerically calculated to be e^l'i/sw), where 
I^SVk) are the eigenfunctions of Hsw- We now com- 
pare these eigenfunctions with the numerically evaluated 
exact eigenfunctions of HuaU ■ 

In order to make the comparison, we calculate the fi- 
delity between the approximate eigenfunctions and the 
exact ones. The fidelity is defined by: 



/=l(* 



approx I ^ exact 



(C5) 



For various system parameters, the fidelity for the first 
few eigenvectors are plotted in Fig. |4l In Figs. 4(a) 
and |4(b)[ we choose the system parameters to be such 
that the qubit frequency is much smaller than the oscil- 
lator frequency. In these quasi-degenerate qubit cases, we 
see that the eigenvectors calculated within the adiabatic 
approximation have better fidelity than the eigenvectors 
evaluated according to the SW transformation. On the 
other hand, as can be seen in Figs. 4(c) and |4(d)] when 
the qubit frequency is more than the oscillator frequency, 
the adiabatic approximation breaks down completely but 
the SW transformation method still holds. 

Thus, in the quasi-degenerate qubit regime, the adi- 
abatic approximation is better than the SW transfor- 
mation in approximating the eigen-structure of the Rabi 
Hamiltonian. In this paper, we only explore the quasi- 
degenerate qubit regime. For this reason, we carry all 
our analysis based exclusively on the adiabatic approxi- 
mation. 
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